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Abstract 

We introduce a one-parameter deformation of the Wishart-Laguerre or chiral ensembles of 
positive definite random matrices with Dyson index j3 = 1,2 and 4. Our generalised model has 
a fat-tailed distribution while preserving the invariancc under orthogonal, unitary or symplcctic 
transformations. The spectral properties arc derived analytically for finite matrix size N x M for 
all three /?, in terms of the orthogonal polynomials of the standard Wishart-Laguerre ensembles. For 
large- N in a certain double scaling limit we obtain a generalised Marccnko-Pastur distribution on 
the macroscopic scale, and a generalised Bessel-law at the hard edge which is shown to be universal. 
Both macroscopic and microscopic correlations exhibit power-law tails, where the microscopic limit 
depends on (3 and the difference M — N . In the limit where our parameter governing the power-law 
goes to infinity we recover the correlations of the Wishart-Laguerre ensembles. To illustrate these 
findings the generalised Marcenko-Pastur distribution is shown to be in very good agreement with 
empirical data from financial covariance matrices. 
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1 Introduction 



Ensembles of matrices with random entries have been introduced in the pioneering works by Wigner 
and Dyson in the early 1950s, in an attempt to provide a statistical description of energy levels 
of heavy nuclei based on few and minimal symmetry requirements. Today's applications cover all 
areas in Physics, and we refer to pQ for a review. However, many years before Wigner and Dyson, 
John Wishart had already introduced a random matrix ensemble in his studies on samples from 
a multivariate population. This ensemble of random matrices can be obtained as follows: take a 
rectangular matrix X of size M x N (M > N), whose entries are independent Gaussian variables 
with mean zero and variance one in the real, complex or quaternion domain, labelled by the Dyson 
index (3 = 1, 2 and 4, respectively. Then form the (non-normalised) covariance matrix of size N x N, 
W = X^X, which is positive definite. By comparing to such matrices one can establish which part 
of the spectrum of empirical covariance matrices carries genuine information and which should be 
discarded as contaminated by pure noise. 

The ensemble of such matrices W is known in Random Matrix Theory (RMT) as Wishart, Laguerre 
or chiral ensemble and will be denoted by WL in the following. It has since appeared in many different 
contexts, such as multivariate statistical data analysis [2], analysis of the capacity of channels with 
multiple antennae and receivers [3], low-energy Quantum Chromodynamics and other gauge theories 
[H [5], Quantum Gravity [6j [7], knowledge networks [Sj, finance [9] and also in statistical physics 
problems, such as a class of (1 + l)-dimensional directed polymer problems [10] . Very recent analytical 
results include statistics of large deviations for the maximum eigenvalue [TTJ and distributions related 
to entangled random pure states |12j . 

Because of its invariance, the WL ensembles can be written in terms of eigenvalues and completely 
solved, both for any finite N and M, as well as in the large- N limit. We only give a brief account on 
what is known, more details follow in the main part later. 

The WL ensemble can be exactly solved through the standard orthogonal polynomial technique 
|13j . For a Gaussian weight functions all spectral correlation functions can be expressed at finite- N in 
terms of orthogonal (j3 = 2) or skew orthogonal (/? = 1,4) Laguerre polynomials, see [13] for a review. 
In the large- N limit one has to distinguish between smooth macroscopic, and oscillatory microscopic 
correlations on the level of mean spacing between eigenvalues. The most known macroscopic result 
in the spectral density of Marcenko and Pastur (MP) [14] . generalising the semi-circle in the Wigner- 
Dyson class (WD). On the microscopic domain much effort has been spent on the so-called hard edge 
at the origin for the positive definite eigenvalues, which is absent in the WD class. It is governed by 
the Bessel-law [151 HI EE E3 H] labelled by (3 and M — N being finite. The robustness or universality 
under polynomial deformations of the Gaussian weight function has also been proven [191 120] . This 
perturbation destroys the independence of the uncorrelated degrees of freedom of X, without changing 
the microscopic correlation functions. Connected macroscopic density correlation functions also remain 
universal under such deformations [6]. 

What is the motivation to further generalise WL ensembles? In the finance and risk management 
domain, the empirical covariance of a set of N assets over a temporal window of size M has been 
under scrutiny for some time [9], and its eigenvalues were shown to be distributed in reasonably good 
agreement with the MP law, as if they were originated by a completely uncorrelated data series. 
However, the same analysis repeated by several groups [211 Ell [23] on different data sets have shown 
that either the part of the spectrum corresponding to extremely low eigenvalues - the most interesting 
for portfolio selections - or the fat tails are not reproduced by this crude approach. This has led to 
the appearance of more sophisticated models [23], [24J, [25] , e.g. the multivariate student distribution 
where the variance of each matrix entry becomes a random variable. Another deformation of the 
WL ensemble has been introduced in [26] where sparse matrices X were considered. This setting has 
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many applications in communication theory and in complex networks (namely in the study of spectral 
properties of adjacency matrices). All these generalisations lead to a deformation of the MP law, and 
thus lie outside the WL class. However, the lack of invariance in such models generically spoils the 
complete solvability for all correlations functions. 

Conversely, a model with power-law tails which is exactly solvable in principle appeared in |25j . 
However, the analysis in was restricted to the macroscopic spectral density, whereas many more 
interesting results, and ultimately a complete solution of the model can be obtained. This goal is 
achieved by exploiting techniques and results introduced in previous works, where the WD class was 
generalised using non-standard entropy maximisation J27J l2H] and super-statistical approaches [29l [30] . 

We will follow these lines and provide a complete solution for all three /3 of a generalised WL model 
with rotational invariance, with an emphasis on the issues of universality and complete integrability 
for all spectral correlations, both macroscopic and microscopic. After a proper rescaling and normali- 
sation, our A-independent correlations depend only on a single parameter a, given by the power-law. 
In the limit of a large parameter a > 1 we recover all standard WL correlations. For completeness, 
we also mention another generalisation of WL which exploits a different direction. In [31] the unitary 
WL were generalised to display critical statistics. 

Our paper is organised as follows. In the next section [2] we define our generalised WL ensembles as 
a one-parameter deformation, including a general non-Gaussian potential V for all three /3 = 1, 2 and 
4. The general solution for finite-A is given applying the method of (skew) orthogonal polynomials 
to an integral transform of the standard WL ensembles. 

In the next section [3] we take the macroscopic large-A limit for the special case of a Gaussian 
potential. Here macroscopic refers to the smooth part of the spectrum, considering correlations on 
a distance large compared to the mean level-spacing. We rederive a generalisation of the semi-circle 
and MP spectral density, see refs. \27\ [28j [29] and [25] respectively. As a new result we compute the 
average position of an eigenvalue and the position of a pseudo edge. The former leads to the correct 
scaling with A and an A-independent generalised spectral density displaying a power-law decay. 

Section U] is devoted to the microscopic large-A limit, where correlations on the order of the mean 
level-spacing are computed. Here we can use the known universal WL results as an input to our 
model. In the case when the difference of the matrix dimension M — A is kept finite the spectrum has 
a hard edge at the origin. In two subsections we compute its microscopic density there, generalising 
the universal Bessel-law as well as the corresponding smallest eigenvalue distribution for all three (3. 

In section [5] the nearest neighbour spacing distribution in the bulk of the spectrum is computed 
using a Wigner surmise at N = 2 for our generalised model. 

Our conclusions and outlook are presented in section [6l including a comparison of the generalised 
MP density to data from financial correlation matrices. In two appendices technical details are col- 
lected. 

2 Definition of the model and finite- N solution for general potential 

The joint probability density of our generalised Wishart-Laguerre ensembles is defined as follows in 
terms of matrix elements 



where X is a matrix of size M x N with real, complex or quaternion real elements for the values 
(3 = 1, 2 or 4, respectively. We define M = N + u for later convenience, where v > may be either 
finite or of order O(N) in the large-A limit. The integration measure c£X is defined by integrating 
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over all independent matrix elements of X with a flat measure. Expectation values of an operator O 
(denoted by (C) 7 ) are defined with respect to V~, in the usual sense. 

The real positive parameter 7 is to be specified below, and we keep an additional variance-like 
parameter n > 0. The so-called potential V is taken to be a polynomial of finite degree d, although 
some of the universal results we inherit from the WL ensembles are known to hold for a much larger 
class of functions. 

The measure f)2. 1 1) is well-defined and integrable only if the following condition holds 

7 d > ^-N{N + u) . (2.2) 

This condition, which is derived in appendix [Aj can be seen when changing to radial coordinates for 
the matrix X^X, see e.g. |32[l28j. In particular, the large- N behaviour of any spectral property cannot 
be taken for fixed 7, and a prescription about the way both quantities should approach infinity needs 
to be given, respecting the inequality eq. (12. 2ft . 

A similar model was introduced earlier generalising the (non-chiral) Wigner-Dyson ensembles for 
a Gaussian potential [271 [29J [28], and a similar interplay between the deformation parameter 7 and 
the matrix size N was observed. 

In the limit of an infinite deformation parameter 

lim TUX] = exp[-n/3TrF(X t X)] = P[X] (2.3) 

7— >oo 

we recover the standard Wishart-Laguerre or chiral ensembles denoted by PfX]. Typically, in these 
ensembles one chooses n ~ N when taking the large- N limit to obtain an A-independent macroscopic 
density. In our case, this limit will be more involved, and we will keep n general for the time being. 

The generalised ensembles can be related to the standard WL through the following integral 
representation: 

(l + z)^ = —rj o dC e"* p- 1 e-t* . (2.4) 
Inserting this into the definition eq. (|2.ip we obtain 



1 r°° 



-^W(XtX) 
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(2.5) 



This relation is the key starting point to solve our generalised model both for finite- and large-A. 
The same trick was used for the generalisation of the Gaussian Wigner-Dyson ensembles introduced 
previously in |27 [ I28 [ |2"5], In fact, a similar technique was employed much earlier in [32] when solving 
the fixed and restricted trace ensembles by writing them as integral transforms of the Wigner-Dyson 
ensembles. 

An advantage of our model over some other generalisations of WL [25[ [23] is its invariance under 
orthogonal, unitary or symplectic transformations. In particular, for (3 = 1 Burda et al. [25] considered 
a very general family of probability distributions of the form: "P/[X] ~ /(TrX T C _1 XA~ 1 ) where C 
and A represent the true correlation and autocorrelation matrices respectively, and / is a non-negative 
and normalised weight function. Only in the special case C = A = 1 invariance is recovered. This 
approach has been modified in [23] to allow for a time-dependent random volatility. 

From eq. (|2.5|) (or eq. (|2.1[) ) we can immediately go to an eigenvalue basis of the positive definite 
matrix X^X, to obtain the following joint probability distribution function (jpdf) 

V ^i=l / i=l j>k 

= f^J Q d£e-H 1 - 1 V(\ 1 ,...,\ N ;Z). (2.6) 
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It is expressed through the jpdf of the standard WL 



N 



V(X\, . . . ,Xn;0 = II A 



|j9(i/+l)-l 



cxp 
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j>k 



(2.7) 



depending on £ through its weight exp 
from the integration over the angular degrees of freedom 

For completeness we also define the corresponding partition function 



. In both jpdf 's we have suppressed the constant 
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j>k 



r(7) 
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which is again an integral over the standard, ^-dependent WL partition function 

N a r „ 1 N 



2(0 = 



o 



\d\i A, 



cxp 



i=i 
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j>k 



(2.9) 



Because of this linear relation between the generalised and standard ensembles we can immediately 
express all fc-point eigenvalue density correlation functions, denoted by R for finite- N, in terms of each 
other. They are defined in the usual way [13] 



Hy(Xl, . . . , Afc) 



1 



(N - k)\ Z. 
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C 5 r( 7 )^ 



dXk+i ■ ■ ■ dXwV^Xi, ... , Xn) 
R(X 1 ,...,X k ;0 , 



(2.10) 
(2.11) 



where the fe-point correlation functions of the standard ensembles depend on £ through the exponent 
in the measure 



R(X 1 ,...,X k ;0 



(N-k)\ 1 



dX k+ i ■ ■ ■ dX N V(Xi, . . . , Ajv; £) . 



(2.12) 



The latter can be solved using the method of (skew) orthogonal polynomials [13], expressing them 
through the determinant of the kernel of the orthogonal polynomials for (3 = 2, or the Pfaffian of the 
matrix kernel of skew orthogonal polynomials for (3 = 1 and 4. We only detail the simpler (3 = 2 case 
here and briefly outline (3 = 1 and 4, where we refer to [13] for more details. 

Let us define monic orthogonal polynomials and their norms for (3 = 2 as follows 



dX X u exp 



P fc (A)fl(A) = h k 5, 



kl ■ 



(2.13) 



Introducing their kernel 



K N (X,fi) = (Xn)2e 



N-l 



h^PkWHOj.) > 



(2.14) 
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and applying the Christoffel-Darboux identity for A ^ fi 

2^ h k PkWPkW = h N-i V— " » ( 2 - 15 ) 

fc=0 " 

we can express all eigenvalue correlations of the standard ensemble through this kernel [13] , 

R(X 1 ,...,X k ;0= det [K N (Xi,Xj)]. (2.16) 

l<i,j<k 

We thus arrive at 

Ry(\i,...,\k) = r^^^W^TiM ■■AK N (X i ,X j )] , (2.17) 

which is the main result of this section. The simplest example is the spectral density i? 7 (A) given by 
the integral over the single kernel Kn(X, A). Notice that it is normalised to N = J °° dXR^(X). 

In order to take N large we only need to know the asymptotic of the polynomials Pjv, take a finite 
determinant of size k of the asymptotic kernel and integrate once with respect to £. In appendix [B] the 
orthogonal polynomials, the corresponding densities and partition functions are worked out in detail 
for finite- N and a Gaussian potential V(X) = X at = 2, given in terms of Laguerre polynomials and 
their norms. 

As pointed out already the same result eq. 112. 17j) holds for = 1 and 4 when replacing the 
determinant by a Pfaffian, Pf[Kjv(Aj, A,-)], where kat is a 2 x 2 matrix kernel. For a Gaussian weight its 
skew orthogonal polynomials are explicitly known as well in terms of Laguerre polynomials [171 [T8] . 

For completeness, we also give the partition function occurring inside the integrand in eq. (|2.17l) 
in terms of the norms hi of orthogonal polynomials 

N-l N-i 

z(o = m n hi = m h$ n r ?~ i > ( 2 - is ) 

i=0 i=0 

or their ratios n = An identical result holds for = 1 and 4 in terms of the skew orthogonal 

norms |13j . 

In the Gaussian case the ratio of partition functions Z(£)/T(7)iJ 7 in eq. ([2.17P can be obtained 
most explicitly at finite- N for all 3 values of by changing to radial coordinates, see appendix lAl for 
a derivation. 



3 Macroscopic large- N limit for a Gaussian potential 

In this section, we will restrict ourselves to the Gaussian potential V(X) = A, deriving a generalisation 
of the MP spectral density of WL for all three 0. It exhibits a non-compact support and power-law 
behaviour for large arguments, and we will distinguish two cases. In the first subsection, we will deal 
with matrices that become asymptotically quadratic with M — N = v = 0(1). In the standard WL 
ensembles, at large-iV the spectral support is an interval on the positive semi-axis, where the origin 
represents a hard edge. In the second subsection, we take the limit 

N . . 

lim — = c , (3.1) 

M,N^oo M 

with c < 1, corresponding to the case M — N = u = 0(N). In the WL ensembles, the resulting MP 
macroscopic spectral density takes support on a positive interval. Conversely, in both cases c = 1 and 
c < 1 our generalised macroscopic density will have support on the full positive real semi-axis. 
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3.1 Generalised semi-circle for c = 1 



The finite- N density for the WL ensembles R(X; £) is well known in terms of Laguerre polynomials, 
see eq. (|B.7p for f3 = 2 in the appendix [Bj This results into an explicit integral representation for the 
spectral density of our generalised ensembles, see eq. ([B.lOp , Despite this result, it is rather difficult 
to extract information about the macroscopic large- N limit from those analytical formulae, both for 
the standard and generalised WL ensembles. 

Hence, we follow here an alternative route, already exploited in [271 128} [29]: we directly insert the 
large- N result into eq. (|2.1ip . The N ^> 1 asymptotic of the WL density is known and given by the 
Marcenko-Pastur law, which is in fact the semi-circle law in squared variables at c = 1 



n 2N 

lim R(X) = -\ — - 1 , with A e (0,2JV/nl . (3.2) 
7V>1 7T V n\ 

It is given here for the Gaussian weight exp[— n{3\] for all three /3, and can be derived easily using 
the Coulomb gas approach and saddle point method (see also [15]). 

In order to obtain an iV-independent macroscopic density we have to rescale the argument of the 
density by the mean eigenvalue position, A — ► {\)x, and divide by N to normalise the density to unity. 
The mean position of an eigenvalue or first moment, (A) 7 , can be computed in the generalised model 
for both finite- AT and finite-^ (see appendix [A]). 



J °° dX A Ay (A) 1 
\ # 7 (A) 
l{N + u) 
2n( 7 - §N(N + v) - 1) 



(3.3) 



We will comment on the existence of this first moment later, and we will also need this equation again 
when we consider c < 1. It correctly reproduces the known result for WL in the limit lim 7 ^ 00 (A) 7 = 
(A) = (N + u)/2n. For the WL case we thus obtain the following known N- and /3-independent 
macroscopic density from eq. (|3.2|) 



p{x) = lim hx)R{x{X)) = ^-J--l, with xg(0,4]. (3.4) 

It is normalised to unity and has mean (x) = 1. 

We can now repeat the same steps for our generalised model, where we will need eq. (I3.2p for the 
weight exp[— £nf3X/^/] (see eq. (|2.7p ). Because of the rescaling with respect to (A) 7 we now have to 
specify the A r -dependence of 7. We keep the following combination fixed, 



a = lim 

7V,7— > 00 



7 - ^N(N + v) - 1 



(3.5) 



with a > finite. This satisfies the constraint (j2.2p . We thus obtain for the generalised macroscopic 
density 

p a (x) = lim -J-(A) 7J R 7 (x(A) 7 ) 

N, 7— >oo iv 



A 1 



hm — A 7 / d£ e — ,/____ - 1 , (3.6 

rt^coN ' J T T(7)2: 7 7vr U <x(A) 7 
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where the integration is restricted to the interval X = (0, ^jrn ] = (0, 4a/ x]. In order to compute the 
integral we still need the following quantity inside the integrand, 

(3.7) 



r (7)^7 r (j - f N(N + u) 



the ratio of partition functions given here for finite- N. This result is derived in appendix |Aj see eq. 
(IA.5j) . Inserting all ingredients into eq. (|3.6p and taking limits with the definition (|3.5p we arrive at 
the following result after changing variables, 



1 /4cA Q+2 f 1 , 

Pa[x) = 2anT(a + l) [t) J dt eXP 

3\ /a„.\<*+2 



4a 
1 

x 



. /I _ i 



= ==^ ^7 r — iFi a + -;a + 3; . (3.8 

4a v / ?r(o + l)r(a + 3) \ x J \ 2' x) y ' 

Here we have introduced the confluent or Kummer hyper geometric function 

iF^a-b-z) = — r(6 . ) TV . /" dt e zt t a - 1 (l-t) b - a ~ 1 . (3.9) 



r(6 - a)r(o) 

Eq. ()3.8|) is the main result of this subsection, the macroscopic spectral density of our generalised 
WL. It has an unbounded support (0, oo), and the density as well as its first moment are normalised 
to unity 

roo roo 

/ dxp a (x) = 1=1 dx x p a (x) . (3.10) 
Jo Jo 

Note that due to this normalisation, the parameter n has completely dropped out. We are left with 
a one-parameter class of densities depicted in fig. [H which approach the WL density for a — > oo as 
discussed below. 

From the expansion for small arguments 1-^1(0; b; z) = 1 + %z + . . . we can immediately read off 
the power law decay of our new density eq. 



hm p a {x) = x-(" +2 > r J { ^f ( ^o, (4a) Q+1 (1 + 0(1/*)) ■ (3-H) 

Because of a > the decay is always faster than quadratic. However, if we drop the requirement for 
the existence of the first moment we can allow for a to take values — 1 < a < while satisfying the 
constraint eq. (|2.2j) . Keeping the same formal rescaling in eq. (13. 6p we arrive at the same result eq. 
(|3.8p now valid for — 1 < a and a 0- Our density can thus describe power laws in between linear 
and quadratic decay as well, see fig. Q3 The same feature could be incorporated in the generalised 
Gaussian Wigner-Dyson model [271 [29| [28] . 

As a check we can take the limit a — > 00 on our final result eq. (|3.8p . This amounts to decoupling 
the 7- and A r -dependence, and thus we expect to recover the MP density at 7 = 00. By taking a 
saddle point approximation, we find that this is indeed the case, liniQ^oo p a (x) = p(x). Hence, the 
density p a (x) is a well-behaved deformation of the MP density for c = 1 . 

We can also derive the behaviour of the density p a {x) close to the origin. Using the large argument 
asymptotic for the confluent hyper geometric function at negative argument, lim^oo \F\{a\ b; — z) = 
z- a T{b)/T(b - a)(l + 0(l/z)), we find that 

Iimp a (x) = x-l r r / a *f* (l + 0(x)) • (3.12) 
x-*o 7rl [a + ljva 



x We use the obvious notation q q+1 = exp[(a + 1) In |a| 
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Therefore all our generalised densities have a square root singularity at the origin, just as the MP 
density. 

To illustrate our findings we first plot eq. (|3.8p in Fig. [T] (left) for different values of a, and 
compare it to the semi-circle density eq. (|3.4p . In order to visualise the square root singularity for all 
a we map the density from the positive to the full real axis by defining 



®<x{y) = \y\p a (y 2 ) , 



(3.13) 



that is to a normalised density on R, dyd a (y) = 1. In this form it equals the deformed semi-circle 
law derived from generalising the Gaussian Wigner-Dyson ensembles in \27\ [2"9"1 [2"5] , where we have 
eliminated all irrelevant parameters. 

The same map eq. (|3.13p takes the MP density eq. (|3.4[) to the semi-circle, as mentioned already 
several times, 

1 



m = s 



ir 



(3.14) 



Both densities are shown in Fig. [T] (right). For a = 14 it already approximates the semi-circle very 
well, whereas for a — ► the generalised density rapidly approximates a ^-function. For negative 
— 1 < a < the height of the maximum in fig. Q]left goes down again, for a = —0.5 the curve is even 
below the semi-circle. 



p a (x) 



My) 



0.8 



. 6 



. 4 



0.2 




Figure 1: The macroscopic generalised density eq. (|3.8|) p a (x) shown on the positive real line R_|_ 
for a = —0.5, 0.1, 0.5 and 14 in light blue, blue, red and green, respectively (left), and its map 
&a(y) = \y\Pa(y 2 ) to the full real line R (right). Note that the MP or semi-circle density given in black 
for comparison has compact support on (0,4] and [—2,2] respectively. 



As a final point of this subsection let us discuss the issue of macroscopic universality of our 
generalised density, eq. (|3.8p . This is a direct consequence of the (non-)universality of the MP or semi 
circular density, due to the linear relationship eq. (12. lip . 

The semi-circle possesses a certain weak universality in the sense that it is the same for all three 
Gaussian ensembles at (3 = 1,2,4, as well as for independently distributed random variables, as was 
shown already by Wigner. Consequently our generalised density is universal in this weak sense, too. 

On the other hand the invariant deformation of the Gaussian potential by a polynomial in the 
definition of our model is clearly non-universal, as the semi-circle becomes a polynomial times one or 
several square root cuts, and we refer to [6] for details at (3 = 2. Remarkably, it was found in [6] for 
(3 = 2 that all macroscopic two- or higher /c-point connected density correlation functions are universal 
under such perturbations V, depending only on a finite number of parameters for any degree d. Does 
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this universality persist in our model? The answer is no, simply by looking at the definition of the 
connected two-point density: 

R™ nn {\,H) = Ry(X, n) - Ry(X)Ry(n) . (3.15) 

We use the same definition for the standard WL density. It no longer relates linearly to the corre- 
sponding connected WL two-point density as we would subtract an integral of the product of two 
1-point densities, instead of the product of two integrated 1-point densities: 

7 Jo r (7)^7 

Jo r (7)^7 Jo r (7)^ 7 

jt r d^e^ C'- 1 -^- R conn {\^;i) . (3.16) 
Jo T Yl) z i 

The universal macroscopic connected two-point function obtained in the large- N limit p conn (x, y; £) in 
[6] will thus mix with the non-universal density p(x; £), and the same feature persists for higher fc-point 
connected correlators. This should not come as a surprise as the same situation was encountered in the 
fixed or restricted trace ensembles \32\ [33] , being an integral transformation of the classical Wigner- 
Dyson ensembles. As observed there |33] our microscopic correlations will remain universal, see sect HI 
In the next section we will study the limit — > c < 1. The corresponding standard WL ensemble 
can be mapped to the so-called generalised Penner model [7] with positive definite matrices. The extra 
determinant from the Jacobian of this change of variables can be written as an extra logarithm in the 
potential V — > V + iVln |A|. For the same reason as given above the universal findings made in [7j for 
the unitary ensemble (3 = 2 do not translate to the macroscopic limit in the next section either. 



3.2 Generalised Marcenko-Pastur law for c < 1 



In this subsection we deal with the limit in which the matrix X remains rectangular, that is both M 
and N = cM become large such that limjv,Af^oo N/M = c < 1. This limit is particularly relevant for 
applications to real data series. 

We will follow the same steps as in the previous subsection. It is known that in the standard WL 
ensembles with weight exp[— n(3\) the corresponding density is given by the following MP law, 



lim R(X) 

AT»1 




N 
2n" 



-X. 



N 
2n" 



-X: 



A 



with A G 



N 
2n' 



-X. 



N 
~2~n~' 



-X: 



(3.17) 



For later convenience we have defined the bounds 

X± = (c~5 ± l) 2 , with < c < 1 . (3.18) 

In the limit c — > 1 we recover from eq. ()3.17|) the semi-circle eq. (j3.2|) from the last section. 

The A r -independent density is again obtained after rescaling by the mean eigenvalue position 

7 N 



(A), 



2nc( 7 - &N 2 - 1 



(3.19) 



where we have used M = N + v = N/c. For large 7 we obtain the quantity (A) = which is the 
average position for the standard WL in our limit c < 1. We thus obtain for the rescaled density MP 
density 



p(x) = lim -(\)R(x(X)) 

N— >oo iv 



1 



2lTCX 



\J [x — cX_){cX + — x) , with x G [cX_ , cX^ 



(3.20) 
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It is normalised to unity with mean (x) = 1. 

For the generalised model we have to insert eq. (|3.17p . now with weight exp[— £— A], into eq. 
(12, lip and rescale with respect to eq. (|3.19j) . As previously we keep fixed 



a 



lim 



7 



2c 



1 



with q > 



as in eq. f)3.5j) . The rescaled generalised density is thus given by 



p a (x) = 



lim -J-(A) 7 i? 7 (a;(A) 7 ) 

N, 7^00 iv 



(3.21) 



(3.22) 



N. 



lim 1(A) 7 / dee-^V 

7— >oo iv 1 



g(0 < 

(7)^7 jnx(\) 




x(A), 



where X = [f£JT_ 



COL V" 



Filling in all definitions and changing variables we finally arrive at the 



following main result of this subsection: 

1 



p a (x) 



2ircaT(a + 1) 



,a+2 f X + 


' ca " 


1 / dt exp 


1 




X - 



t a ^(t-X.){X + -t) 



(3.23) 



Our density in normalised to unity and has first moment (x) = 1. This result was derived previously 
(modulo different notations) for (3 = 1 in [25], using different methods. The integral in eq. p. 231) 
can be computed in principle in terms of a confluent hypergeometric series in two variables (see [34J, 
formulae 3.385 and 9.261(1)), but the integral form is more convenient for numerical evaluations and 
an asymptotic analysis. As a first check we recover eq. (|3.8p in the limit c — > 1. Furthermore, one can 
show that the following limit holds, linic^oo p a {x) = p(x), recovering the MP density eq. (|3.20p . Thus 
the large- iV limit and the large-7 limit are again well behaved. We have also checked in appendix 151 
that the convergence with N towards the density eq. (|3.8p is very fast, see fig. [U in fact faster than 
for WL. Our remark from the previous subsection allowing — 1 < a < applies here too. 
It is easy to see how the density decays for large x S> 1 



lim p a (x) 



>+2). 



(ca) 



Q + l 



27rr(a + 1 

with the same power law as for c = 1. The constant C is given by 

'3 _ X. 



C(l + 0{l/x)) , 



c 



X. 



-a; 3; 



X_ 



The asymptotic for small values of x is less obvious to obtain, and we find 

V , 



lim p a (x) 



x- a - l l 2 exp 



_ca x _ 



(3.24) 



(3.25) 



(3.26) 



where the constant V = X*(X + - X_) l / 2 {ca) a ~ 1 / 2 /\QT(a + 1). 

While the MP density of the standard WL ensemble eq. (|3.20p has compact support between 
[cX-,cX+], our generalised density is non- vanishing on the entire real positive axis, even for c < 1. 
To the left of the edge in MP, x < cX_, our density decreases but remains non-zero. Below a 
certain point that we will call pseudo edge, X—, our density becomes exponentially suppressed. From 
the asymptotic (|3.26|) it is possible to give an estimate for X_, below which the density becomes 
negligible. The reasoning goes as follows: writing the asymptotic (|3.26p as 



p a (x) ~ exp 



ca 



'a + 1/2) log x + — X. 

x 



(3.27) 
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Figure 2: The macroscopic generalised density p a {x) eq. (|3.23p for a = 3 and c = 0.3 (blue), compared 
to the MP distribution p(x) (|3,20p (black). In the inset, the behaviour close to the origin is shown. 
The red dashed line corresponds to the pseudo edge X- of our generalised MP density (see main text 
for details). 

the exponential damping conventionally begins at the point X- where 

coi 

{a + 1/2) log X. + — X_ « 1 (3.28) 

A.— 

For the case depicted in fig. [2] (c = 0.3 and a = 3), the above estimate reads X_ « 0.055..., in 
reasonable agreement with the inset. 

4 Universal microscopic large- N limit for a general potential V 

In this section we consider a different large- iV limit, the microscopic limit, which takes us to the scale 
of the mean level spacing and thus to the distribution of individual eigenvalues. Our findings will be 
universal for a general non-Gaussian polynomial potential V(A) for all three j3 = 1, 2, 4, inheriting the 
corresponding universality from the standard ensembles. 

We will only consider the case c = 1 and the so-called hard edge here, deriving a generalised 
Bessel-law for the microscopic densities and their first eigenvalue distributions. For c < 1 the local 
distribution at the inner (and outer) soft edge of the standard WL ensembles follows the Tracy- 
Widom-law. Although it would be very interesting to derive the corresponding generalisation we have 
not managed so far, and leave this task for future investigation. 

In the first subsection the microscopic densities are derived while the second subsection is devoted 
to the first eigenvalue distributions. The matching of the two is illustrated in many pictures throughout 
this section, being an important consistency check. 
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4.1 Generalised universal Bessel-law 



Let us first recall the definition of the microscopic limit in the standard WL ensembles, resulting into 
the universal Bessel-law. We will discuss in detail the case (3 = 2 and then only quote the results for 
= M. 

For simplicity consider the Gaussian case first. Starting from the weight exp[— n[3\] we first scale 
out the mean eigenvalue A — > x{\), just as in eq. (|3.4p for the macroscopic limit. On top of that we 
make a further rescaling by the mean level spacing 4N 2 x = y, keeping y fixed. We therefore define 
the microscopic limit as 

where R(\) is the spectral density for finite- iV in one of the three WL ensembles. The result now 
depends on (3 and v as indicated through the indices, in contrast to the semi-circle law. We can then 
apply the following asymptotic 

to the orthogonal Laguerre polynomials in the finite- ./V density, see eq. (|B.7p for (3 = 2. We obtain 

= l£dtJ u (Vtyf = l(MVy) 2 - J„-i(Vv)Ju+i (Vy)) , (4.3) 

after replacing the sum by an integral! with variable t = k/N . The Bessel density is plotted for different 
values of v in fig. 2] together with the first eigenvalue from the next subsection, after changing to the 
conventional squared variables y — > y 2 (see eq. f|3. 13j) ) . 

^u\y) = ^{Uy) 2 - Ju~i{y)Jv + i{y)) • (4.4) 

This result is universal [19] being valid for any potential V with spectral support including the origin. 
We only have to rescale in eq. (|4.ip by the macroscopic density in terms of the squared variables 7r$(0) 
for a general potential V, instead of the Gaussian macroscopic density eq. (|3.14p where 7r$(0) = 1. In 
other words all orthogonal polynomials eq. (|2.13p tend to Bessel- J functions in the microscopic limit. 

We can now repeat the above analysis for our generalised microscopic density. The scaling of 7 
with N keeping a fixed is kept throughout this entire section. In the previous section we found that 
the density diverges exactly like the standard density as an inverse square root, see eq. (|3.12p . For that 
reason the microscopic rescaling is the same, without changing powers of N. However, the constant 
in front of the macroscopic density at the origin is not the standard Gaussian result, 1/V = i?(0), but 
given by eq. (I3T2|) . <& a (Q) = b/ir with 

r(a + ^ (4.5) 



T(a + l)\/a 

We therefore define the microscopic limit as 

1 ,u n ( 



2 Instead of eq. (|B.7j) we could have applied the Christoffel-Darboux identity eq. (|B.9[l . leading to the same result. 
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Because of the universality we just stated we can restrict ourselves to do the computation for the 
orthogonal polynomials with Gaussian weight exp[— 2n^X/j] in eq. f)2. 13|) . Taking the microscopic 
limit eq. (|4.6p and inserting eq. (14. 3h we obtain 

This is the first main result of this subsection given here in terms of squared values. The integral 
could be expressed in terms of generalised hypergeometric functions, but for plots this representation 
is preferable. Note that in our the calculation we have inserted the ratio of partition functions eq. 
(|3.7p for the Gaussian models. This quantity is again universal as in the large- N limit not only the 
polynomials themselves, but also their norms become universal. Starting from eq. (|2.18p the standard 
WL partition functions is given as follows 



lim (]ag[Z(£)]-]ag[Nlhg]) = lim N 



N^oo 



N- 



N-l 

£ 

i=0 



N 



logN 



di(l-t)log[r(t)] (4.8) 



The ^-dependent ratios of the norms determined by the so-called string or recursion equation at finite- 
N have a universal limit r(t) [19]. Re-exponentiating and inserting this universal result eq. (|4,8p into 
eq. (|2.8p . the generalised partition function Z^ also becomes universal in the large- N limit, and thus 
the ratio eq. (|3.7p as well. The constant factor that we have subtracted on the left hand side of eq. 
(|4.8p cancels out when taking the ratio. 




(2) 

Figure 3: Varying a at (5 = 2: the generalised microscopic density da' v (y) eq. (|4.7j) (blue) and its first 
eigenvalue (green) at a = 0.1 (left), a = 2 (middle), and a = 20 (right), vs the corresponding WL 

/o\ _____ 

Bessel density 'do (y) eq. (|4.4p (black) and its first eigenvalue (red). 





(2) ______ 

Figure 4: Varying v at j3 = 2: the generalised microscopic density $ a 'u{y) eq. (|4.7j) (blue) and its first 

(2) ______ 

eigenvalue (green) at a = 0.1 vs the corresponding WL Bessel density (y) eq. (|4,4p (black) and its 
first eigenvalue (red): v = (left), v = 1 (middle) and f = 2 (right). It is clearly visible that even the 
first eigenvalue of the generalised model has fat tails. 
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The new microscopic density eq. (|4.7p generalising the Bessel-law eq. (I4.4p is shown in figs. Eland 
d] for various values of a and v. The plots include the corresponding first eigenvalues to be derived 
later. As a check we can take the limit a — > oo to analytically reobtain eq. (|4.4j) from eq. (14. TJ) - This 
is illustrated in fig. [3] where we observe that the convergence is rather slow. We have checked that for 
a ~ 0(150) the first three maxima become indistinguishable. 

The procedure for = 1 and 4 is the same as above and so we can be more concise. The 
corresponding microscopic densities of the WL ensembles are again universal |20j. Computed initially 
in [TTJ, HH] the obtained expressions can be simplified. They can be expressed through the = 2 
density eq. (|4.4p plus extra terms as shown for (3 = 1 [35] and = 4 [55] , 




WHv) = 0{v) + \ju(\y\) [i- I dtJ v (t) , (4.9) 



€\y) = 4i{^y)-\j2u{2y) dtj 2v {t) 



(4.10) 



The generalised densities immediately follow. Because of the linear relationship they are also expressed 
through the generalised (3 = 2 density eq. (|4.7j) : 



poo I t -i /I7/I \ / r\y\~\/i/ a b 2 \ 

€Uv) = # ( &(y) + r(a + 1) y o # ^S.- Jv {m y fi/^ (1 - J du v (t)){A.ii) 
tfgUv) = *2W(2y) - r^T+i) 7 ^ e ^ r V ^2 ^(f^) i ^(t)(4.i2) 



For a given f the inner integral over the single Bessel- J function can be performed analytically. It is 
given terms of Bessel- J functions for odd values of u, e.g. Jq dtJ\{t) = Jq{v), and additional Struve 
functions for even v. 

The generalised microscopic densities i?^i(y) are compared to the standard ones below in fig. [5] 
for (3 = 1, and in fig. [B]for (3 = 4. 

Higher order correlation functions can be computed along the same lines by inserting the asymp- 
totic Bessel kernels into eq. (|2.17p . and we only quote the simplest final result for (3 = 2: 

k 



x det 

l<i,j<k 



y% - Vj 



(4.13) 



The corresponding results for (3 = 1, 4 are given in terms of a Pfaffian of a matrix kernel [To], and for 
a discussion of a relation between the three universal kernels we refer to [20] . 

A feature we observe for all three (3 is that for a < 0(1) the oscillations of the Bessel density are 
completely smoothed out, apart from the first peak. A similar feature was observed in a generalisation 
of the unitary WL ensemble for critical statistics [31]. However, no power law tails seem to be present 
in such a model where only the generalised microscopic density and number variance were computed. 

It is known that for standard WL the maxima of the Bessel density correspond to the location of 
individual eigenvalues [37], as we will see in the next subsection. On the other hand the microscopic 
density of the WL ensembles in the bulk is completely flat, equalling - in our normalisation. We may 
thus suspect that in the generalised model the bulk is approached much faster than in the standard 
WL, where localised maxima persist to i/ > 10. We therefore focus mainly on the first eigenvalue 
distribution in the generalised model which is the subject of the next subsection. 
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4.2 Generalised universal first eigenvalue distribution at the hard edge 

The probability that the interval (0, s] is empty of eigenvalues is defined as follows, 



EM) 



1 



7 Js 



dXi • • • dAjv^ > 7(Ai, . . . , Ajv) 

e .7-1 ^(0 



r( 7 )z 7 

where the gap probability of the WL ensembles is defined as 

1 



2(0 



d\ 1 ---d\ N V(\i,...,\ N ;£) . 



(4.14) 
(4.15) 

(4.16) 



Both quantities are normalised to unity at s = and vanish at s 
eigenvalue p(s) simply follows by differentiation. 



pM) 



d_ 

ds 



EM) , 



oo. The distribution of the first 



(4.17) 



and likewise for WL. In WL the gap probability E(s) and the first eigenvalue distribution p(s) are 
explicitly known and universal in the microscopic large- N limit for all v at (3 = 2, for odd values of 
v and at (3 = 1, and for v = at (3 = 4. This has been shown by various authors independently 
|16[ [38, 39, 37]. In some cases only finite- N results are know in terms of a hypergeometric function of 
a matrix valued argument |40[ I41j. from which limits are difficult to extract. 

Although p(s) follows from E(s), the most compact universal formulas are known directly for p(s) 
for all three (3 [37] . There, also the second and higher eigenvalue distributions are given, which we 
will not consider here. 

For pedagogical reasons we start once more with an explicit calculation for (3 = 2 and the Gaussian 
ensemble. At v = the pre-exponential factor is absent and we have for WL with weight exp[— 2nA] 



1 f c 



d\\ ■ ■ ■ d\w exp 



N 



-2n^ Ai 



i=l 



N 

j>k 



A fc p 



exp[— 2nNs] 



(4.18) 



shifting all integration variables by s and using the invariance of the Vandermonde determinant. This 
result is exact for any TV and identical to the properly rescaled large- N result when keeping Ns fixed. 
The generalised ensemble follows easily, by inserting this result into eq. (|4.15p 



EM) 



exp 



2n£Ns 

7 



1 + 



2nNs\ 
7 J 



-( 7 -7V 2 ) 



(4.19) 



which is also exact for finite and infinite N. This very fact implies that we have full control of the 
large- N limit. Because of the fat tails in our distributions it was not a priori clear if the scales at 
the hard edge and in the bulk would mix. This implicit assumption of a separation of scales in the 
previous subsection is thus fully justified. 

The microscopic limit can be taken following eq. (14. ID for WL: 



lhn E(-^r 



(A) 



(4.20) 



where we explicitly indicate the dependence on (3 and as in the previous subsection. As a result we 
obtain for (3 = 2 and v = 



7(2) 

'u=0 



(y) 



exp 



1 

"4 y 



(4.21) 
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and for the corresponding generalised gap probability 



S. 



(2) 



(y) 



lim E, 



y 



4N 2 b 2 



1 + 



\ub 2 



-(a+l) 



(4.22) 



The first eigenvalue distribution can be compared to the microscopic densities $ a (y) m squared vari- 
ablcs: 



p?2 (y) 



dy to [V ] 



1, 

^\y\ ex P 



1 2 



for WL, and for the generalised ensemble 



p£l=o(y) 



dy a > o[y ] 



{a + l) 
2ab 2 



1 + 



Aab 2 



-(o+2) 



These distributions are all normalised to unity, 



i 



(4.23) 



(4.24) 



(4.25) 



The restriction to a Gaussian potential in the discussion above can be lifted, as the first eigenvalue 
distributions for all v are universal, including the ratio in partition functions that we have inserted 
again. Eqs. f|4.23[> and (|4.24j) are compared to the corresponding densities eqs. f)4.4j) and (|4.7p in fig. 
[3] for various values of a. 

Next we give the first eigenvalue distribution for general v. Here we directly use the most compact 
universal expression [37] for p^\y) in WL, without making the detour over £^ 2 \y) [39] . 



P?Hy) 



i, 

2^1 ex P 



1 s 



det [Ii- j+2 {\y\ 

l<i,3<v 



(4.26) 



In addition to the exponential in eq. (j4.23[) it contains a determinant of finite size v x v over the 
modified Bessel-/ function, which is absent at v = 0. Knowing that the properly rescaled microscopic 
gap probability is a function of the form = (y/ £y / ab 2 ) , see eq. (14.71) . we obtain for the 

generalised first eigenvalues distribution in terms of squared eigenvalues 



1 r 00 f) 



1 



r(a + i) y 



* 2aP 



cxp 



4a& 2 



det 



(4.27) 



For v = \ containing only one Bessel-I the integral can be performed and is given in terms of a 
hypergeometric function 



p { £i(y) 



2ab 2 T(a + 

r(q + 3) 
|y|r(a + l) 



/•oo 



^ 2 



4a& 2 



4a 



-(a+3) 



4a \ -i 

1+ y 2 " 



(4.28) 



It is shown in fig. [5] (middle), together with v = 2 (right). For increasing v however, the integral 
representation eq. (j4.27H is more convenient. 
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Figure 5: Varying v for [3 = 1: the generalised microscopic density (y) eq. (liTTj) (blue) with its 
first eigenvalue (green) at a = 0.1 vs the corresponding WL Bessel density $^}\y) eq. (|4.9p (black) 
and its first eigenvalue (red) at v = (left), v = 1 (middle), and v = 3 (right). 



We now turn to (3 = 1. Here the first eigenvalue distribution of the WL ensembles is only known 
explicitly for odd values of v in the large- N limit [37], given by a Pfaffian with indices running over 
half integers: 



p$\y) = const. \y\W 2 



cxp 



1 



Pf_ f + i<i, m -i [(i - j)I i+j+ 3(\y\)] ■ (4.29) 



The constant in front is determined by the normalisation to unity and can be computed case by case. 
An exception is v = where the distribution was calculated in [16] 



pg\y) = -l(2 + |y|)exp 



1 



\v\__y_ 

2 8 



For v = 1 and 3 we have from eq. (14.291) 



pi\y) 



^|yl ex P 



exp 



1 



h{\y\) ■ 



(4.30) 

(4.31) 
(4.32) 



The corresponding generalised formula for general odd v thus reads 

1 



p£l(v) 



const, 



(3-u)/2 



exp 



8ab 2 



x Pf -f+i<*.i<f-i 



(«- VIA*) 



(4.33) 



up to the normalisation constant. For the simplest examples v = 0, 1 and 3 displayed in fig. [5] we have 

(i), N i r A , -u<*+\( Ay\ [I\ \ \y\ fl Zy 2 



Pa,o(y) 

p2s(v) 



4^r( c 

(a + l) 
4a6 2 



8afr 2 



, (4.34) 

(4.35) 
(4.36) 
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While the density is modified only rather mildly compared to WL, the first eigenvalue changes 
considerably. We have checked that the curves converges to WL for large a, where the convergence to 
the density is much faster than for (3 = 2. The case v = in fig. [5] left is the only example where the 
microscopic densities do not vanish at x = 0. The fact that they both have the same limit \ can also 
be seen analytically, exploiting that J u (0) = 5 Vt o (see eqs. (|4.9p and (|4.1ip ). As mentioned above for 
v = 2k with k S N+ the first eigenvalue is not available to date. 

Finally we turn to j3 = 4. In principle the result is known in the WL ensemble, 



Po\v) = const. \y\ 



3 

"2 exp 



1 2 

-2 V 



z z/2{{\y\hv 



(4.37) 



Here Z 3 / 2 ({y}2i/+i) is the large- N matrix model partition function at topological charge 3/2 with 2u+l 
degenerate masses at value y, and we refer to [37j for a more detailed discussion of these objects. This 
partition function is generally known explicitly only for an even number of masses, except at v = 0. 
There we have 



P ( o\y) 



1 , . .3 

- v 2ir \y\ 2 exp 



1 



h/2(y) = \y\ (cosh(y) 



1 

y 





1 2 


) exp 


~2 y _ 



(4.38) 



Thus the generalised distribution depicted in fig. [6] is given by 

£ dC ^ r+1 (|,| cosh (fv^) - ^sinh 



P { o\y) 



aT(c 





[ iy 2 ~ 




exp 


2ab 2 






(4.39) 




Figure 6: Varying v for (5 = 4: the generalised microscopic density da^{y) eq. (|4TT2l) (blue) at a = 0.1 
vs the corresponding WL Bessel density -d^\y) eq. (|4.10p (black) at v = (left), v = 1 (middle), and 
v = 2 (right). For u = we also display the respective first eigenvalue for the generalised (green) and 
WL ensemble (red). 

We note that the oscillations get smoothed out as already observed at (5 = 2. The convergence in 
a to the WL quantity is again slow, as for (3 = 2. 



5 The Wigner surmise in the bulk 

In this section we study the spacing distribution between eigenvalues in the bulk in our generalised 
model. In contrast to the previous sections we do not take iV to be large, but use the N = 2 results, 
following the original idea of Wigner. 

In the Wigner-Dyson ensembles with a Gaussian potential, the spacing distribution has the simple 
form 

pW(s) = As^ 2 , (5.1) 
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known as Wigner's surmise. The known constants A, B follow from normalisation and can be found 
e.g. in pp. For the WL case the corresponding expression can be computed from the jpdf, where we 
introduce v = ^(y + 1) — 1: 

VW{ 8 ) = Cs^ +1 / 2 K 1/2+i? (nPs) (5.2) 
where K^x) is a modified Bessel function and the constant C is given by: 

C = U-y^^in^-^-^V 0-^) r (l + v + f ) ) 1 (5-3) 

Only for D = 0, one recovers the WL Wigner's surmise, using i^ 1/ / 2 (^) = \/tt/2x e~ x (apart from 
squared variables in the exponent there). For general v one easily gets an exponential decay for large 
spacing ~ s /3+v e -n/3s ^ f rom the asymptotic expression of the Bessel function, lim^^oo K fl (x)^ir/2x e~ x . 
However, we expect that only the v = expression will lead to a good approximation of the infinite- iV 
case. 

Conversely, for the generalised model the spacing distribution has a different expression: 

pfM = ^ + |)" + " 7 ,fi (-P, p + l; -P + 7 ; i - jX.) (5.4) 

where 2-^1(0, fr; c; z) is a hypergeometric function. The constant C 7 can be computed as 

2 1 -^(P + 1,-217-1 + 7 )r(7)r(l + g/2) (nf5\ 2 + 2 ^^ 

7 r(i? + i)r(i + p)T(v + 1 + /9/2)r( 7 - 2 - /? - 2*) V 7 / 1 J 

and B(a,(3) is a Beta function. The large-s decay is now a pure power law ~ where we have 

defined 

w = i-P-2u-2 > 0, (5.6) 

required to be positive for convergence^]. 

For v = 0, the spacing distribution takes a much simpler form: 

'~< [s) ~ 7 r03)r( 7 - p - 2) V 7 ) ' ( j 

It agrees with the corresponding quantity in the generalised WD ensembles found in [29J. In this 
particular case, we can compute the mean level spacing explicitly: 



00 



«.» 7 = I .*$«(.)* _(-L)-l±^, (5.8) 
which converges to the WL Wigner surmise value for 7 — > 00: 

((,»= r sT (P)( S ) ds =±±p. (5.9) 

JO n P 

After defining the rescaled quantities having mean spacing 1, 

Tf (x) = ((s)),V^ (« S » 7 *) (5.10) 

V^\x) = ((s»P (/3) (((s»x) (5.11) 



3 This condition equals eq. (|2.2|l derived for the partition function at N = 2. 
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we can compare the curves for all three /3 at v = (and n = 1), in fig. The power-law tail 
modification compared to the standard WL spacing distribution is evident in the plots. Because in 
the large- N limit 7 scales with N, we keep the combination w in eq. (15. 6ft fixed to be able to compare 
to spacing distributions at large- N. 

Both in this and the previous section the power-law tail of the macroscopic density is seen to 
persist on the microscopic level of the mean level-spacing. It would be very interesting to confirm this 
on real data sets. 




Figure 7: Comparison between Vj (x) (blue, green) and VW( X ) (red), for = 1,2,4 (from left to 
right). The 7 value for the blue curves is chosen is chosen in such a way that the combination w 
is kept constant to the value 2. The green curves have value 7 = 12, 12, 25 from left to right, and 
correctly approach the limiting WL curve. 



6 Conclusions and outlook 

We have introduced a generalisation of all three ensembles of random matrices called Wishart-Laguerre 
or chiral ensembles, replacing the exponential of a non-Gaussian potential by a fat-tailed distribution 
with parameter 7. In the limit 7 — > 00 we can recover the exponential weight and thus the standard 
ensembles. This modification lead to the appearance of correlations with a power-law, governed by 
a single parameter. Such a behaviour is found in many systems in nature, e.g. in the wide area of 
complex networks. WL ensembles are often used in comparing to eigenvalues from covariance matrices 
of real data sets. To illustrate the potential of our generalisation we show a comparison to financial 
data in fig. The eigenvalues of the covariance matrix from time series of stock data clearly show a 
power-law behaviour. These are well described by the generalised Marcenko-Pastur density, refining 
previous comparisons to the standard Marcenko-Pastur law. 

The solution of our generalised models relied heavily on the possibility of writing them as an 
integral transform of the standard WL ensembles. This generalisation is thus in the spirit of super- 
statistics where other models have been constructed already. The virtue of our model is its invariance 
for all three symmetry classes, allowing to go to an eigenvalues basis and to study universality. We 
could show that the generalised macroscopic density which was known in the Gaussian case is only 
weakly universal. 

In contrast all microscopic densities are universal under any invariant deformations by polynomial 
potentials. This macroscopic/microscopic dichotomy should not come as a surprise, being observed 
previously for the restricted trace ensembles. 

We exploited the linear relation to standard WL to solve our model exactly at finite-iV for any 
polynomial potential, using the formalism of orthogonal polynomials therein. In the subsequent large- 
N double scaling limit, where 7 is scaled with N, we derived all density correlations in the macroscopic 
limit for quadratic and rectangular matrices, and in the microscopic limit at the hard edge for all 



21 



three values of (5. Here we have mainly focused on the spectral density itself and the first eigenvalue 
distribution. The general formalism for higher density or higher individual eigenvalue correlation 
functions was provided and is straightforward to use if such quantities will be needed. 

While the hard edge was solved exhaustively, persisting in our model for asymptotically quadratic 
matrices, we only provided a Wigner surmise in the bulk. Here, more detailed properties of correlation 
functions could be investigated, including a possible generalisation of the Tracy- Widom distribution 
at the soft edges. This is left for future investigations. 




0.5 1 1.5 2 2.5 3 3.5 



Figure 8: Comparison between the rescaled eigenvalue distribution from financial data [42] and the 
macroscopic density p a (x) for the generalised model eq. (|3.23p , in red dots and solid blue respectively. 
The best fit gives a value of a ~ 0.95, which corresponds to a power-law decay as p a (x) ~ x -2 ' 95 . 
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A Partition function and first moment of the Gaussian models 

The purpose of this appendix is threefold. First, we derive the condition under which our generalised 
model defined in eq. (|2.1[) is convergent. Second, we compute both the generalised and standard 
Gaussian partition functions in order to determine their ^-dependent ratio needed in the computation 
of all eigenvalue correlation functions. Third, we compute the first moment as a function of N, u, (3 
and 7 which is needed for the rescaling of the eigenvalues in the large- N limit. 

All three steps will be performed by changing variables from independent matrix elements or its 
eigenvalues to radial coordinates, following [32]. 
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The generalised partition function reads in terms of eigenvalues 

N N , N 
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Here we have changed to radial coordinates of the iV-component vector of the eigenvalues (Ai, . . . , Ajv)> 
and da,N denotes the angular integration over the ^-dimensional unit sphere Q(N). The eigenvectors 
e« = Aj/r of norm unity span Q(N) and no longer depend on the radius. Collecting all powers of r in 
the numerator and comparing to the leading power of the denominator 07 at large r>l, the integral 
only converges if the following inequality holds: 







N(N + v) -I- < -1 , 



(A.2) 



which is exactly eq. (|2.2p . 

Next we compute the partition function where for the rest of this appendix we restrict 

ourselves to the Gaussian potential V(A) = A. The same steps can be taken for a purely monic 
potential V(A) = X d as well. 

In principle we could repeat the same calculation in terms of eigenvalues as above, but it will be 
more instructive to start directly from the matrix elements: 
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(A.3) 



Here we have used radial coordinates for the (5N{N + v) component vector of all independent matrix 
elements Xjj, with squared norm r 2 = TrX^X. We don't need to compute the angular integral 
explicitly as it cancels out below. 

If we insert the result eq. (|A.3P into the relation (|2.8|) we immediately obtain 
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Combining the last two equations we arrive at 
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the ratio of the two Gaussian partition functions at finite values of N, v and 7 valid for all three 0. 
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In the last step of this appendix we compute the mean eigenvalue position in the Gaussian model. 
It can be either defined through the spectral density, see eq. (j3.3|) . or in terms of the first moment, 
where we start with the standard WL ensembles: 
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(A.6) 



Note that the /3-dependence has cancelled out. We can immediately use this result to compute the 
same quantity for the generalised Gaussian model, 
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^Tr (XtX) 
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(A.7) 



This result is used in the rescaling of both large- N limits. As a check it reduces to eq. (|A.6j) in the 
limit 7 — ► oo, with a weight exp[— n(3r 2 ]. 



B Explicit j3 = 2-solution for all /c-point densities at finite N and 7 

In this appendix we present all details for the solution of the generalized WL ensemble with unitary 
invariance (3 = 2 and Gaussian potential V(A) = A. In this case the orthogonal polynomials of the 
WL ensemble are know to be Laguerre, allowing for an explicit solution at finite ./V and finite 7. 
The orthogonal polynomials in eq. (12.130 read for the weight function exp[— £^A] 

aw = (-)**■ (^)'^(M A ), (B „ 

with norms 

/ \ 2k+v+l 
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Here the Laguerre polynomials are defined as usual 



L%{z) = E(-iy(^ > with L l^' = - L £l&- (B-3) 

j=o ^ 3/3- 



Using eq. (12.180 we can immediately read off the WL partition function from the norms, 



N Nv+N 2 N-l 
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As the next step we can compute the partition function eq. (|2.8p given by 
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This leads to the following ratio needed for example inside the relation (|2.1ip 
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It confirms independently part of the result from the previous appendix, eq. (IA.5D for (3 = 2. 

The spectral density for finite A follows by inserting this ratio as well as the standard Laguerre 
density at finite- A, 
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With the help of eq. flEjj) we can derive and simplify the Christoffel-Darboux identity for Laguerre 
polynomials of equal arguments 
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We thus arrive at our final result for the generalised density at finite A and 7: 
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This single integral over an exponential times polynomials can be performed explicitly, at the expense 
of a double sum. Since this equivalent result is not very illuminating or useful for the asymptotic we 
do not display it here. 

Proceeding along the same lines as above we can write down the general result for the /c-point 
density correlations functions as they follow from eq. (I2.17P 
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In order to compare the finite-A result ([B~TU]) with the macroscopic 7- and A-independent density 
i? Q (x) eq. (13. 13f) . we adopt the following procedure: 
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1. We rescale i? 7 (A) with mean value (A) 7 eq. (|3.3|) and normalise to 1: py(x) = N 1 (A) 7 i? 7 ((A) 7 x). 

2. Next, we express 7 as a function of a and A, 7 = a + N(N + v) + 1, and pass to squared 
variables: #a(x) = |x|p 7 (x 2 ). 

3. Then, we compare $ a (x) and ^^(x) for v = in Fig. [9j The agreement is already very good for 
N = 4, apart from the region close to the origin. 
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Figure 9: The macroscopic generalised semi-circle density $ a (x) eq. (|3.13|) for a = 1.02 and 14 (green 
and dashed red), compared with the finite N = 4 result (blue and dash-dotted orange). 
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